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ABSTRACT 

Topological properties of the vacuum were investigated in finite-temperature lattice 
QCD without and with dynamical quarks using a local operator of topological 
charge density and a variant of cooling. Below the deconfinement temperature 
evidence for flux tube formation was observed in the distribution of topological 
charge density around a static QQ pair. With dynamical quarks, the flux tube 
breaks when the distance between the static quark and antiquark increases. 



1. Introduction 

Many phenomenological models of strong interactions motivated by QCD share a 
set of simple assumptions: confinement of quarks, antiquarks and gluons; nontriviality 
of the QCD ground state; existence of a flux tube between Q and Q with constant 
energy density per unit length; existence of a phase transition temperature Tc above 
which confinement and fiux tubes disappear. 

Lattice QCD is a powerful tool to test the above assumptions nonperturbatively 
from first principles. The existence of fiux tubes has been proven in numerous lattice 
studies both indirectly, by showing that the QQ potential rises linearly for large QQ 
separations, and directly, by measuring correlations of plaquettes with the Wilson 
loop or a pair of Polyakov lines. ^'^ 

Colour confinement is believed to be due to the structure of the QCD vacuum, 
which in turn should be a complicated superposition of quark and gauge fields with 
nontrivial topology. Then it is natural to expect that the existence of the fiux tube 
between a QQ pair will also be refiected in quantities characterizing topological prop- 
erties of vacuum configurations containing a QQ pair. The aim of the present contri- 
bution is to show that the existence and basic properties of fiux tubes can be inferred 
from distributions of the topological charge density around a static QQ pair. 

2. Topological charge on a lattice and cooling 

The topology of a continuum gauge field configuration is characterized globally 
by the (integer) value of the topological charge Q: 
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Q = ld'x q{x) = -A_ld'^ e^^^^ f;^{x) f;,{x), (1) 

where q{x) is the topological charge density. A natural choice for defining Q on a 
lattice is to replace the integral in Eq. 1 by a sum over all lattice sites, Ql = J2x ^l{x)^ 
with some lattice operator of topological charge density qL{x) which in the naive 
continuum limit converges to q{x). Here we will restrict ourselves to results which we 
obtained using the plaquette definition"^: 

1 ±^ 

(1l{x) = - ^432^2 Yl ^A^^p^ Tr [U^^x) UpAx)], (2) 

with U^,y being a plaquette with (/i,z/)-orientation.* Employing this definition, a prob- 
lem immediately arises: topological charges of lattice configurations are non-integer 
and in general rather small. The origin of this problem was clarified by Campostrini 
et al}: lattice and continuum versions of the theory are different renormalized quan- 
tum field theories which differ from one another by finite renormalization factors, e. g. 
(1l{x) = a'^ Z{fi) q{x) -\- 0{a^). Z{fi) is rather small for the range of fi values used in 
Monte Carlo simulations. 

A way to get rid of renormalization constants while preserving physical informa- 
tion contained in lattice configurations is the cooling method. ^'^ One assumes that 
small local changes of a configuration do not change its global topological properties. 
We applied the variant of cooling suggested by Hoek et al7 

3. Results 

Our results come from three runs, each of them performed on an 8"^ X 4 lattice 
using the Metropolis algorithm. In the first two runs we evaluated path integrals 
in pure SU(3) gauge theory with standard Wilson action for fi = 5.6 (confinement 
phase) and 5.8 (deconfinement phase). The third run was done for QCD at /3 = 5.2 
with SU(3) Wilson action for gluons and 3 fiavours of Kogut-Susskind quarks of equal 
mass ma = 0.1 (confinement phase again). We made lOOk iterations in each run and 
measured our observables after every 50th sweep. Each recorded configuration was 
subjected to 50 cooling steps. 

In all simulations, cooling quickly removed the renormalization factor Z{fi) and 
the measured values of Ql converged (close) to integer values. We then observed sub- 
stantial differences between lattice configurations below and above the deconfinement 
phase transition. The Ql = sector dominates in both phases, but in the confine- 
ment phase almost 70% of configurations have non-zero charge and the distribution 
of charges is approximately Gaussian (in both pure QCD and full QCD), while in the 
deconfinement phase Ql 7^ appears only very rarely. 

The most interesting result comes from measuring (iv(O)iv^(J) qlir)), the correla- 
tion function which is related to (the square of) the topological charge density around 
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Fig. 1. The square of the topological charge density around a static quark-antiquark pair, 
{L(0)L^(d) q^(r))/({L(0)L^(d)){q^)), at distances 3 and 4, after 5 cooling steps. 



a static QQ pair with separation d [L denotes the usual Polyakov hne). In both phases 
we see a suppression of the topological charge density around static charges, consis- 
tent with phenomenological models. However, there is an essential difference between 
the phases which is illustrated in Fig. 1. In the deconfinement phase the suppression 
exists only in the vicinity of static sources (lb,e), while in the confinement phase 
it occurs in the whole fiux tube between Q and Q, both without (la,d) and with 
dynamical quarks (Ic). In the pure gauge case the fiux tube is most clearly visible at 
the largest separation, J = 4 (Id). In full QCD for J = 4 (If) we see an indication of 
fiux tube breaking due to creation of a virtual qq pair. A similar effect was observed 
earlier also in colour field distributions around a static QQ system.® 
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